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Abstract - 

THa libration of a gravitationally oriented rigid satellite on 
an e llip tic orbit is studied by the method of canonical transformation. 
The possibility of the center of mass of a satellite traveling on 
an elliptic orbit is examined. The transformations are performed in 
a mathematical vay according to an operational method of perturbation 
conserving Hamiltonian formulism. The orbit eccentricity effect 
upon the amplitude and the frequency shift of the libration are ■ •* ' 
exhibited. The ultimate nature of satellite libration on an elliptic 
orbit is revealed.. £ 
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Nomenclature - 
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X, Y, Z 


Pr in cipal, central moments of inertia of the satellite 
Semi-major axis. « . 

Eccentricity. 

Hamiltonian. ' 

Integral constant, 
lament of inertia about f 
Inertia ~ensor. 

Unit vectors along principal axes. 

Gravitational parameter of the center of force., 
Langrangian. 

-Action variable. 

Angle variable. 

Mass of satellite . 

Unit vector along 

Momentum conjugated to q. .. ■ 

Angle defined by equation (15). 

Initial value of q. 

Radius from the center of force to the center of mass 
Potential energy. 

Kinetic energy. y 

Eccentricity perturbation functions. i 
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6 True anomaly.* 


<? 

a 


Angle of libration . 

Constant Js U~ ^'X' J 





Introduction - 


This paper is devoted to studying the planar li brat ion of a 
gravitationally oriented rigid satellite on an elliptic orbit. The 
equations describing the notion of the satellite are derived using 
the Lagrange method. The coupling of the orbital and librational 
motion does not permit the satellites center of mass traveling on 
an elliptic orbit. By introduction of canonical variables the 
librational motion assumes a form particularly suited for the 
Von Zeipel technique. The mathematical treatment and the physical 
consideration in the present approach produce a solution vhich reveals 
the ultimate nature of satellite libration. The application of the 
canonical transformation method Cl) to this type of problem is apparently 
new in the literature on satellite dynamics (2-lh). Based on this 
analysis,- a proposed damp ing method for satellite libration will be 
discussed in a later article . 



Equations of Motion - 


A satellite of mass m with principal axes ji, and k at its 
center of mass is launched into an elliptic orbit. The longitudinal 
axis i. and the transverse axis £ lie in the orbital plane and the 

i * 

axis k is normal to the orbital plane. The inertia, tensor at the 
center of mass is 


( 1 ) 

where A, B and C are principal, central moments of inertia of the 
satellite about axes i, and k respectively. 

Let (y, 6) be the polar coordinates of the satellite’s center of 
mass and cp which is in the orbital plane be the angle of inclination 
of the longitudinal axis _i to the line joining the center of force to 
the center of mass. \ \ 

The kinetic energy of the satellite is 


/A 0 0\ 
J = o 3 0 

\ O O cl 
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By MacCullsgh’s formula, the potential energy Is 



Km 
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+ 0(-fc) 
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in which K is the gravitational parameter of the center cf force. 
Thus for motion about the earth K = ©4e where G is the constant of 
gravitation and Me is the mass of earth. I » n • ^ • n and n is a 
unit vector along y* Accordingly, . 


Z = A Crffr BSSf 
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So second order in 9 the Lagrangian of the motion is 
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The three equations which govern the three degrees of freedom S 


of the satellite motion are 




-where h is an integral constant, a is the semi-major axis and e the 
eccentricity of the orbit, equation (8) becomes 



By substituting equation Cll) in equation ( 7 ), the only solution to 

these equations is <p *» 0. This is expected physically because if 

<P ^ 0 the center of mass of the satellite is moving in a non-Newtonian 

C 

force field which is not central. However, if — is much less than 

m 

unity, it is acceptable to examine the libratlonal and orbital motion 
separately. • • ’} . 

Substituting equation ( 9 ) into equation (7) the result is 
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Let the independent variable be transformed from the time to the 


polar coordinate, true anomaly 0. (tee fjnds that 


•/: 'I 
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Accordingly, equation (12) becomes 


cfr_ zes±e_ di. > /-f fM)d, 

dfr JteCsi& 4& T 3/t €(^hc. <Z 


Uk) • 


/mate 




/t-ecae 
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After introducing 


■mr 


equation (l^) reduces to 


v •!' . 
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dP . rjf/-eV( it) r **** U _ 
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where p 


= 


is the momentum conjugate to q. 


Equation (l&) can be expressed In the form 


tszy =r £2? ^ 


vhereln 2 = 3(l-e 2 )(— ~— ) > 0. For circular orbit, i.e-, e = 0, the 

w 

librations will evidently be purely harmonic with constant frequency 
B-n^ the’ amplitude of libration will be constant over time depend in g bn 
the initial conditions. With the existence of the eccentricity perturbation 
the solution of equation ( 17 ) is either almost periodic or diverging or 

. t ' ... 

decay ing librational motion. A solution which can reveal the ultimate 
nature of the librational motion governed by equation (17) is desired. • 
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^ Canonical Transformation. 

■ ' * * • t i • ‘ 

The canonical equations corresponding to equation (17) are 
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(19) 


H is the Hamil tonian of the. unperturbed, system and X is the perturbing 

function. In thi s expression, perturbation terms of order e 3 ana higher 

« 

i 

are omitted. 

The solution for circular orbit subject to the transformed initial 

* 

conditions q( 0 ) = and p^ Q j ® 0 can be expressed in periodic forms 


^ = ^ CoS £2.6 


P^-_&n.s;«ne 


( 21 ) 


Define l Wle and L s § 60 . tliat equations (20) and (21) are 


transformed to . 


f v 


JzL S2 


(25) 
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Since a transformation to new action and angle variables has been 
given in the form 

=/i'oL. D-tf'u, -) = SZL 

the question is, under what circumstances will the transformation be 
canonical, i.e., under what circumstances the new equations will be 



This transformation of variables can be done in the following 
•- manner. Multiply equations (l8) and (19) by - *|^ respectively and 
add. The sum is \ . 



■> i j/i •zrT ~ > <» -v '/ r \ 17 r « -a- \ , 
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It equations (lS) emd 19 ) are multiplied by f|, - -g ami similar 


operations are carried out the result is 
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equations ( 26 ) and ( 27 ) "6600116 



Hie Lagrangian bracket [X,L] in equations (28) and (29) 6 qua l s unity 
"because the Hamiltonian is unchanged. Therefore, comparing equations 
(28) ,and (29) with equations (24) and (25), one obtains 



Differentiating equations ( 22 ) and (23) the results are 


'e3- _P_ and pT e5>s*J 

d! -S2L ~~Z2p- 

£* 



By substituting these results equations (30) and (31) became 
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General Solution 


( 55 ) 


. i • 

Thus the solution for the librational motion about an elliptic orbit 


is described by 
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where L said l are given by 




The physical interpretation of these results is as follows. The 
amplitude and frequency of satellite librations on an elliptic orbit can 
be determined from the librations in a circular orbit by considering that 
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-Vg g — -r K-g' frequency^ arc varying cue oQ rhe e c cenur a c iuy pa r uur s at i on 


Eh.e variation in amplitude is inherent in equation (3*0. Ehe variation 
in frequency is given by equation (35)* It should be noted that the last 
tern in equation (34) represents secular eccentricity effect if 8 is equal 
to l- Accordingly, the amplitude will build up without bound in case 
q s i. For 0^0, the results show that the amplitude to have no secular 
terms. It does not increase or decrease vithout bound but is represented 
by a constant plus trigonometric terms having frequencies which are 
fractions of the fundamental frequency of the system. • 

By excluding the periodic perturbations equations (34) and (35) reduce 


to 
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Taking these values in order to obtain a first approximation to the 
solution one obtains 




c?V 



Q 




(Therefore the first approximation of the librations is 




Sweetie 


( 36 ) 





For D * 1 , equation ( 17 ) is reduced to a fern with constant coefficients 
and the integral solution increases without hound. According to the 
physical sense of the moments of inertia and in view of the instability 
conditions this is referred to as the resonance case. 

Equation (56) describes the librational motion of artificial 
satellites. The ultimate nature of the libration is concisely revealed. 

The eccentricity effect upon the amplitude and the eccentricity frequency 
shift are exhibited. It can be shown that the solution represented by 
equation (36) meets the 'limiting condition imposed. For circular orbit, 
i.e., e s 0, equation (36) reduces to a solution which abounds in the 
literature of satellite dynamics . Based on this analysis, a proposed 
method which damps the librations and prevents librations from becoming 
arbitrarily large will be discussed in a later article. 
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